Superconducting and charge-density wave instabilities in ultrasmall-radius carbon 

nanotubes 
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We perform a detailed analysis of the band structure, phonon dispersion, and electron-phonon 
coupling of three types of small-radius carbon nanotubes (CNTs): (5,0), (6,0), and (5,5) with diam- 
eters 3.9, 4.7, and 6.8 A respectively. The large curvature of the (5,0) CNTs makes them metallic 
with a large density of states at the Fermi energy. The density of states is also strongly enhanced 
for the (6,0) CNTs compared to the results obtained from the zone-folding method. For the (5,5) 
CNTs the electron-phonon interaction is dominated by the in-plane optical phonons, while for the 
ultrasmall (5,0) and (6,0) CNTs the main coupling is to the out-of-plane optical phonon modes. 
We calculate electron-phonon interaction strengths for all three types of CNTs and analyze pos- 
sible instabilities toward superconducting and charge- density wave phases. For the smallest (5,0) 
nanotube, in the mean-field approximation and neglecting Coulomb interactions, we find that the 
charge-density wave transition temperature greatly exceeds the superconducting one. When we in- 
clude a realistic model of the Coulomb interaction we find that the charge-density wave is suppressed 
to very low temperatures, making superconductivity dominant with the mean-field transition tem- 
perature around one K. For the (6,0) nanotube the charge-density wave dominates even with the 
inclusion of Coulomb interactions and we find the mean-field transition temperature to be around 
five Kelvin. We find that the larger radius (5,5) nanotube is stable against superconducting and 
charge-density wave orders at all realistic temperatures. 



The discovery of carbon nanotubes Q] has lead to a re- 
newed interest in the study of Id electron systems. The 
difference between semiconducting and metallic large- 
radius nanotubes may be typically understood by quan- 
tizing the circumferential momentum of the electronic 
states in a single graphene sheet (see, for instance, Q). 
Less conventional properties of nanotubes include Lut- 
tinger liquid behavior of metallic nanotubes found in 
tunneling experiments (see [j| and references therein), 
Coulomb effects 0, Kondo physics [3, and intrinsic su- 
perconductivity observed in ropes [6j and small-radius 
nanotubes in a zeolite matrix |2J. The main focus has 
traditionally been on the effects of the Coulomb inter- 
action between electrons. However, the electron-phonon 
interaction has also received considerable attention both 
experimentally 0, and theoretically 

[nl[H[i3. Most 

theoretical analyses of electron-phonon interactions in 
nanotubes assume the phonon frequencies to be the same 
as in a graphene sheet and calculate the electron-phonon 
coupling strength from a simplified tight-binding model 
for the 7r orbitals of the C atoms. Such an approach, how- 
ever, may not be suited for ultrasmall nanotubes (such as 
the ones in 0), for which the curvature of the nanotube 
leads to strong hybridization of the a and 7r orbitals, 
which results in a qualitatively different band structure 
[l3| . phonon spectrum, and electron-phonon interactions. 

In this Letter we present detailed analysis of three 
representative small-radius nanotubes, the (5,0), (6,0), 
and (5,5), and discuss possible CDW and superconduct- 
ing instabilities of these systems. The (5,0) nanotube 
is the likely candidate structure for the superconducting 
behavior seen in Q with transition temperature mea- 
sured around 15 K. The radii of the CNTs in this exper- 



iment has been determined to within 0.2 A by Raman 
spectroscopy 0]. We demonstrate that even though 
standard electronic structure approaches for calculating 
phonon frequencies, such as the frozen-phonon approx- 
imation (FPA), run into divergences intrinsic to mean- 
field calculations in Id, they can be analyzed from the 
point of view of the random-phase approximation (RPA) 
for the electron-phonon system and parameters of the ef- 
fective Frolich Hamiltonian can be extracted. The main 
results that we obtain are: (i) The strongest electron- 
phonon coupling for the (5,0) and (6,0) CNTs is to the 
out-of-plane phonon modes. This is in contrast to the 
larger radius CNTs which have strongest coupling to 
the in-plane phonon modes as predicted by the nearest- 
neighbor tight- binding model |2(. (ii) Even when the 
residual Coulomb interaction between electrons is ne- 
glected, the larger radius (5,5) CNT remains stable down 
to extremely low temperatures. For the smaller radii 
(5,0) and (6,0) CNTs, when the residual Coulomb in- 
teraction is neglected, the CDW instability was found 
to be dominant over superconductivity for both types of 
nanotubes. (iii) We include the residual Coulomb inter- 
action between electrons following a model developed in 
Ref. [i^. For the (6,0) CNT we find that the CDW tran- 
sition is essentially unaffected by including the Coulomb 
interaction and we obtain Tcdw=5 K. By contrast, for 
the (5,0) CNT both CDW and SC are suppressed but now 
superconductivity becomes dominant with Tsc around 1 
Kelvin. 



The interaction between conduction electrons and vi- 
brations of a crystal lattice is commonly described by 
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using the Frolich Hamiltonian 

kra qfJ- 

+ 9kTk>T>»cl Ta c k , T , a (a qiJ , +al qfl ). (1) 

krk' t' afj, 

Here cJ. r(T creates an electron with quasimomentum k in 
band r with spin c, a qfl creates a phonon with lattice 
momentum q and polarization /z, and q = k — k' modulo 
a reciprocal lattice vector. The energies of electron quasi- 
particles and phonons in the absence of electron- phonon 
coupling are given by e kr and f2^ respectively and the 
electron-phonon vertex is given by gkrk'r'a^- 

To compute the quasiparticle energies £fc r of the 
representative CNTs, we use the NRL tight-binding 
method 0] which has been tested extensively and pro- 
vides accurate results on a variety of materials. After 
fully relaxing the structures with respect to the atomic 
coordinates, the band structure is calculated. We find 
the band structure predicted by zone-folding to agree 
very well with the calculated band structure of the larger 
radius (5,5) CNT. However, for the smaller radius (5,0) 
and (6,0) CNTs there was found to be qualitative dif- 
ferences as shown in Fig. ^ While zone-folding argu- 
ments predict the (5,0) nanotube to be insulating, the 
band structure clearly exhibits metallic behavior. The 
inner band (with the smaller kp) is doubly degenerate 
while the outer band (with the larger kp) is nondegen- 



erate where we have the exact relation 2kp 



fcf. The 



failure of the zone-folding procedure is due to the strong 
curvature effects, which lead to considerable band shifts 
in small-radius nanotubes, as discussed originally in 0] 
for (6,0) nanotubes based on density-functional theory 
calculations. As a result of these band shifts, for the 
(5,0) nanotubes we have a system close to a Van Hove 
singularity which has a density of states of 0.16 states / 
cV / carbon atom. For comparison, the density of states 
for the (5,5) nanotube is only 0.028 states / eV / car- 
bon atom, so we expect that instabilities of the electron- 
phonon systems for the ultrasmall nanotubes are strongly 
enhanced compared to larger radius nanotubes. 

The electron-phonon vertex given in Eq. can be 
expressed as 
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where M krk , T ^ = ^(ip k r\(V^ ~ F )|^ T '>. Here, V qfl 
is the crystal potential under the presence of a phonon 
specified by the ionic displacements SR n i — ue %qRn e qpi (i) 
and Vo is the crystal potential at equilibrium. We cal- 
culate the magnitude of these matrix elements for the 
coupling between electrons on the Fermi surface to all 
phonon modes. For the (5,0) and (6,0) CNTs we find that 
the strongest coupling are to out-of-plane modes. More 





FIG. 1: Left: The band structure of a (5,0) nanotube where 
we set Ef = 0. Right: The phonon modes that have 
the strongest electron-phonon coupling. Shown is a cross- 
sectional slice of the nanotube containing 10 atoms. Modes 
(a), (b), (c), and (d) are out-of-plane optical, out-of-plane 
breathing, in-plane optical, and in-plane acoustic modes re- 
spectively. The X's and O's denote vectors in and out of the 
page. 



specifically, the strongest overall coupling was found to be 
to the out-of-plane optical mode followed by the breath- 
ing mode which are shown in Fig. ^ This is in contrast 
to the larger radius (5,5) CNT which has dominant cou- 
pling coming from an in-plane optical mode. We point 
out that in general, the phonon modes of CNTs cannot 
be classified as in-plane or out-of-plane 0]. However, the 
modes that have the strongest electron-phonon interac- 
tions for the CNTs we study still allow such characteriza- 
tion (see Fig. HI) . Moreover, CNT dynamical matrix cal- 
culations [l7| show that the eigenvectors of these modes 
are essentially the same as in the graphene sheet. 

We point out that, in general, the phonon mode eigen- 
vectors will be influenced by curvature effects, and will 
differ from the graphene results. We checked the relevant 
modes by using the CNT dynamical matrix, and found 
that they agree well with the zone- folding results ^3 ■ 

Now we move on to the issue of how to calculate the 
bare phonon frequencies f2° in the Frolich Hamiltonian. 
In the standard FPA the frequencies are given by 
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(3) 



where u is the amplitude of the displacement, and A_E cos 
and A_E s i n are the energy differences per unit cell between 
the distorted and equilibrium lattice structures where the 
distortion corresponds to the real and imaginary parts of 
5H, n i = ue tqRn e qfJl (i) respectively. The phonon dispersion 
curve for the out-of-plane optical mode obtained from the 
FPA for the (5,0) CNT is shown in Fig.[U Unit cell sizes 
of up to 400 atoms in increments of 20 (which is the 
number of atoms in the smallest possible unit cell) were 
used in the FPA, requiring the phonon wave vectors to 
be commensurate with the chosen supercell. This mode 
shows giant Kohn anomalies at 2k p for the inner and 
outer bands. 

It is important to realize that the divergence of Q qfl 
shown in Fig. [21 does not imply the divergence of f2° in 
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FIG. 2: The phonon dispersion for the out-of-plane optical 
mode in the (5,0) CNT showing logarithmic divergences at 
2kF- We show the fit to the out-of-plane optical mode for 
which the bare frequencies are extracted as discussed in the 
text. The X's denote values of q at which the calculated FPA 
frequencies were imaginary. 



the Frolich Hamiltonian Eq. In the FPA, the phonon 
frequencies are calculated after the electron-phonon in- 
teraction in Eq. Q have been included, which gives 
anomalous softening at 2Uf due to the well-known Peierls 
instability of electron-phonon systems in Id. We have de- 
veloped a technique to extract the bare phonon frequency 
from the numerically computed f2 q/J using a connec- 
tion between the frozen-phonon approximation and the 
RPA treatment of the giant Kohn anomaly which is ar- 
ticulated in Ref. Briefly, the dressed phonon frequen- 
cies f2 gA1 and the bare phonon frequencies fij 3 , satisfy the 
equation 



(fV) 2 = (rC) 2 + 2^ A(«) 



(4) 



where n M (c7) is the phonon self-energy. Using the calcu- 
lated FPA values for fl qf j,, and the calculated electron- 
phonon coupling values to determine IL_i (q) (in the 
random-phase approximation), we can extract the bare 
phonon frequencies by assuming that they have the form 
fijL = A + Bq + Cq 2 and performing least-squares fitting 
where A, B, and C are adjustable parameters. Using this 
method, we have calculated the bare phonon frequencies 
of the representative nanotubes, thereby constructing the 
effective Frolich Hamiltonians for these systems. 

Now we consider possible instabilities of the electron- 
phonon system. Using the RPA analysis on the Frolich 
Hamiltonian, we find the CDW transition temperature is 
given by 



T, 



CDW 



= 4epe" 



-1/AcDW 



(5) 



Here, Acdw is the dimensionless CDW coupling constant 
given by A C dw = Iffg^P^C )/^^ where ^(0) is the den- 
sity of states per spin at the Fermi energy of the band 





(5,0) 


(6,0) 


(5,5) 


Acdw 


0.26 


0.12 


0.024 


Asc 


0.57 


0.12 


0.031 


|UCDW 


0.24 


0.0019 


0.013 


Msc 


0.19 


0.16 


0.093 


^CDW (K) 


160 


5.0 


7 x 10~ 14 


Ts°c (K) 


64 


0.071 


1 x 10~ 12 


Tcdw (K) 


1 x 10~ 18 


5.0 


2 x 10" 43 


T SC (K) 


1.1 







TABLE I: The CDW and SC coupling constants and transi- 
tion temperatures for the CNTs studied. Asc does not include 
Coulomb screening or the temperature dependent softening 
of phonons discussed in the text. T"c DW gc was computed 
without the residual Coulomb interaction while the residual 
Coulomb interaction is included in the calculation of Tcdw.sc- 
Tsc includes the temperature dependent renormalization of 
A SC (see Eq. JTSJ). 



that is undergoing the transition. In table [I] we sum- 
marize the results for the CDW coupling constants and 
transition temperatures for the nanotubes we study. The 
leading CDW instability for the (5,0) and (6,0) CNTs 
is from the out-of-plane optical phonon mode which is 
shown in Fig. ^ while the leading CDW instability for 
the (5,5) CNT is to an in-plane optical mode. 

To study superconductivity, we use the McMillan for- 
mula 



rp («) 

TsC= L20 eXP 



1.04(1 + A SC ) 



Asc - MscC 1 + 0-62A sc ) 



(6) 



where the dimensionless superconducting coupling con- 
stant is given by 



^=4 ^ 5(e kT )5{e k , T ,)\g kTklT> \ 2 -^-. (7) 

fcTft'T'A, W 



In the above equations, v(0) is the density of states at 
the Fermi energy per spin, (O) = 1400 K [19j is the log- 
arithmically averaged phonon frequency, and /ig C is the 
Coulomb pseudopotential which we will set to zero for 
the time being. The results for the superconducting cou- 
pling constants and transition temperatures is also sum- 
marized in tabled for the nanotubes we study. 

From this analysis, we find that both the CDW and 
SC instabilities for the (5,5) nanotube occur below ex- 
perimentally realizable temperatures. This leads one to 
expect that conventional CNTs of larger radius also be 
stable down to very low temperatures. For the (5,0) and 
(6,0) CNTs, the CDW instability was found to be domi- 
nant which occurs from coupling to the out-of-plane op- 
tical phonon mode. 

Now we consider the consequences of introducing the 
residual Coulomb interaction. We point out that since 
the charge density is not evaluated self-consistently in 
the tight binding method we use, the Hartree term which 
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opposes the formation of the charge-density wave is omit- 
ted in our frozen-phonon calculation of frequencies. This 
term essentially gives the Coulomb energy cost of form- 
ing a non-uniform charge density. Including the Coulomb 
interactions properly should lessen the divergences found 
at 2fcp in the phonon spectra. Introducing the residual 
Coulomb interaction changes our effective Hamiltonian 
to 



The residual Coulomb interaction comes up in more 
subtle ways when considering the superconducting insta- 
bility. By properly dressing the electron-phonon vertices 
as well as the phonon propagator in Migdal's expression 
for the electronic self-energy we find that the renormal- 
ized contribution to the superconducting coupling con- 
stant is given by 



- P h 



where 



- E 

kk' qrr' aa' 



V iA A 



-qr a 



Ck' 



'Cki 



(8) 



(9) 



and 7i -ph is given by Eq. Q). For the Coulomb inter- 
action between conduction electrons, we take the form 
used by Egger et al. in Rcf. 15 which, in position space, 
is given by 



V(r- 



e 2 /n 



^( 2 ;-x') 2 +(2i?sin(^))' 



(10) 

Here, the y-direction is chosen to be along the perimeter 
of the CNT and x measures the distance along the CNT 
axis. A measure of the spatial extent of the p z electrons 
perpendicular to the CNT is given by a z ss 1.6 A and R 
is the CNT radius. For the dielectric constant due to the 
bound electrons, we will take the value k « 2 predicted 
by the model of Ref. [H - 

Including the Coulomb interaction in our RPA anal- 
ysis of the CDW instability, we find that the transition 
temperature is modified to 



T, 



CDW 



= 4epe 



— 1/ (AcDW— MCDw) 



(11) 



where /^cdw = v T (0)V q= 2k F and 2/ T (0) is the density of 
states per spin at the Fermi energy for the band that is 
undergoing the instability. We thus find that the effective 
coupling is directly reduced by including the Coulomb 
interaction. 

With our model for the Coulomb interaction, we find 
for the (5,0) CNT that /icdw = 0.24 which dramati- 
cally reduces T CD w to around KT 18 K. Thus the CDW 
instability for the (5,0) CNT is essentially removed by 
taking into account the residual Coulomb interaction be- 
tween conduction electrons. The case is somewhat dif- 
ferent for the (6,0) CNT, however. For such metallic zig- 
zag nanotubes, the wave functions at — k and k close to 
the Fermi energy correspond to symmetric and antisym- 
metric combinations of atomic orbitals in the graphene 
sheet [l2] . Orthogonality of these wave functions within 
the unit cell of the CNT leads to the significantly smaller 
Mcdw = 0.0019 which essentially does not affect the 
CDW transition temperature. 




K ' \(l-VqXo(q)r 



where A^ is the unrenormalized contribution for a spe- 
cific process of wave vector q coupling points on the 
Fermi surface by phonon mode [i and xo (<?) — (/fe+g — 
fk)/{£k+q —£k) [Uj- The first factor describes renormal- 
ization of the electron-phonon vertex by Coulomb inter- 
action and tends to decrease Tsc while the second factor 
corresponds to phonon softening due to the giant Kohn 
anomaly in Id which tends to increase Tsc- The tem- 
perature dependence described by Eq. (|T2")) is similar to 
the two parameter RG analysis presented in Ref. [2l| . In 
addition to the renormalization of Asc we also have the 
direct repulsion between conduction electrons which is 
taken into account through the Coulomb pseudopotential 
as shown in the McMillan formula Analysis based on 
Eq. (|lf)|> we find that the inclusion of /Xg C eliminates su- 
perconductivity in the (5,5) and (6,0) CNTs [HI- For the 
(5,0) CNT, the main contribution to the renormalization 
of Asc comes the 2kp coupling to the out-of-plane optical 
mode discussed earlier. Taking into account the temper- 
ature dependence in xo(l)i using Eqns. (J5J) and l|12fl we 
find a self-consistent solution of Tsc = 1-1 K. Thus we see 
that inclusion of the Coulomb interactions makes super- 
conductivity dominant over the CDW in ultrasmall (5,0) 
CNTs. We point out that our estimates give a mean- 
field value of Tsc- Below this temperature we expect a 
gradual decrease of resistivity, which may be described 
by the Lang er-Ambegaokar-McCumber-Halperin formal- 
ism [22, l2-'"il| . Discrepancy between this value of Tsc and 
the one observed experimentally of 15 K pfl should not 
be a reason for concern because of the exponential de- 
pendence of the superconducting transition temperature 
on the Coulomb interaction strength and the known dif- 
ficulty in calculating the latter accurately. For instance, 
if we replace our estimated value of /ig C = 0.19 by the 
commonly used (jZq — 0.10, we find a self-consistent so- 
lution for the superconducting transition temperature for 
the (5,0) CNT of T sc = 13 K. 
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